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We make a comparative analysis of some recent one-dimensional Langevin models of the acceler- 
ation of a Lagrangian fluid particle in developed turbulent flow. The class of models characterized 
by random intensities of noises (RIN models) provides a fit to the recent experimental data on the 
acceleration statistics. We review the model by Laval, Dubrulle, and Nazarenko (LDN) formulated 
in terms of temporal velocity derivative in the rapid distortion theory approach, and propose its 
extension due to the RIN framework. The fit of the contribution to fourth order moment of the 
acceleration is found to be better than in the other stochastic models. We study the acceleration 
probability density function conditional on velocity fluctuations implied by the RIN approach to 
LDN type model. The shapes of the conditional distributions and the conditional acceleration vari- 
ance have been found in a good agreement with the recent experimental data by Mordant, Crawford, 
and Bodenschatz (2003). 
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I. INTRODUCTION 

Tsallis statistics p] inspired formalism 0, 0, 01 was 
recently used by C. Beck l| to describe Lagrangian 
statistical properties of developed turbulence; see also 0, 
H, H . In recent papers 0, 0, 0] we have made some 
refinements of this approach. The probability density 
function of a component of the Lagrangian acceleration 
of infinitesimal fluid particle in the developed turbulent 
flow is found due to the equation, 

/>oo 

P(a) = / alp P(a\(3)f(f3), (1) 
Jo 

where P(a\(3) is a conditional probability density function 
associated to a surrogate dynamical equation, the one- 
dimensional Langevin equation for the acceleration a, 

d t a = jF(a) + aL(t). (2) 

Here, d t denotes time derivative, F(a) is the determinis- 
tic drift force, 7 is the drift coefficient, a 2 measures in- 
tensity of the noise, a strength of the additive stochastic 
force, L(t) is Langevin source, a delta-correlated Gaus- 
sian white noise with zero mean, 

(£(*)> = 0, (L(t)L(t')) = 28(t-t% (3) 

where the averaging is made over ensemble realizations. 

For constant parameters 7 and tr, this usual Langevin 
model ensures that the stochastic process a(t) defined by 
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Eq. (J2J) is Markovian. The probability density function, 
P(a\j3), of the acceleration at fixed /3, 

/3 = 7/cr 2 , (4) 

can be found as a stationary solution of the corresponding 
Fokker-Planck equation, 

d t P(a, t) = 0„[-7f(a) + a 2 d a ]P(a, t), (5) 

where d a = d/da. This equation can be derived from 
the Langevin equation J5J using the noise © either in 
Stratonovich or Ito interpretations. Particularly, for a 
linear drift force, F(a) — —a, the stationary probability 
density function, d t P(a,t) = 0, is of a Gaussian form, 

P(a\f3) = C(f3)exp[-pa 2 /2], (6) 

where C(f3) = y/ (3/2% is a normalization constant, a S 
[—00, 00]. The function f((3) entering Eq. Q is a proba- 
bility density function arising from the assumption that 
(3 is a random parameter with prescribed external statis- 
tics. 

While it is evident that the three-dimensional Navier- 
Stokes equation with a delta-correlated Gaussian white 
random forcing belongs to a class of non-linear stochastic 
dynamical equations for the velocity field with which one 
can associate some generalized Fokker-Planck equations, 
it is a theoretical challenge to make a link between the 
Navier-Stokes equation and surrogate one-dimensional 
Langevin models for acceleration such as J2| • This model 
is, of course, far from being a full model of the essential 
Lagrangian dynamics of fluid in the developed turbulence 
regime. 

Review and critical analysis of the applications of var- 
ious recent nonextensive statistics based models to the 
turbulence have been made by Gotoh and Kraichnan |13| . 
An emphasis was made that some models lack justifica- 
tion of a fit from turbulence dynamics although being 
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able to reproduce experimental data to more or less ac- 
curacy. A deductive support from the three-dimensional 
Navier-Stokes equation was stressed to be essential for 
the fitting procedure to be considered meaningful. 

In contrast to the fluid particle velocity, the fluid par- 
ticle acceleration, 

a t = ^ = d t Vi + v k d k Vi, (7) 

which incorporates the Eulerian local acceleration and 
nonlinear advection term, can be measured easier by 
using the Lagrangian framework while in the Eulerian 
framework (fixed probe) this requires measurements of 
the velocity Vi, and temporal and spatial velocity deriva- 
tives, d t Vi and dkVi, where du = d/dx k denotes spatial 
derivative in the Cartesian laboratory frame of reference; 
i, k = 1,2,3. In the Lagrangian framework, the Navier- 
Stokes equation can be written as 

a t = -p^ 1 d l p + vd\vi + fi, (8) 

where p is constant fluid density, p is pressure, v is kine- 
matic viscosity, Vi = dtXi is velocity, and fi is forcing. 
Here, Xi = Xi(xok,t) is the particle coordinate viewed 
as a function of the initial value Xi(0) = XQi and time t 
so that the measurement of time series Xi(t) of some in- 
dividual particle by using a fine finite-difference scheme 
allows one to evaluate its acceleration as a function of 
time by using the Lagrangian relation, 

a,i = d 2 t Xi. (9) 

With the initial data points XQi (Lagrangian coordinates) 
running over all the fluid particles one gets a Lagrangian 
description of the fluid flow. Direct analytical evaluation 
of the acceleration from Eq. JHJ is out of reach at present 
so that one is led to estimate it in some fashion. 

The model J5J belongs to a class of stochastic models of 
Lagrangian turbulence and deals with an evolution of the 
acceleration in time which in accord to the Navier-Stokes 
equation is driven by time derivative of the r.h.s. of 
Eq. (JHJ. This type of modelling corresponds to the well- 
known universality (Kolmogorov 1941, Heisenberg 1948, 
Yaglom 1949) in statistically homogeneous and isotropic 
developed turbulence which is expected to occur in the 
inertial range only statistically. Accordingly, the velocity 
and acceleration become random, and one is interested 
in their probability density functions, or multipoint cor- 
relation functions. This is in an agreement with the ob- 
served temporally irregular character of the velocity and 
acceleration of a tracer particle in high- Reynolds-number 
turbulent flows. By the universality, statistics of the ve- 
locity and statistics of the acceleration do not depend on 
statistics of the forcing and chosen initial data. In this 
paper we are interested in statistics of one of the accel- 
eration components, a, so that we model its evolution in 
time. 

In a physical context, an essential fluid particle dynam- 
ics in the developed turbulent flow is described here in 



terms of a generalized Brownian like motion, a stochastic 
particle approach, taking the particle acceleration © as 
the dynamical variable. Such models are generally based 
upon a hierarchy of characteristic time scales in the sys- 
tem and naturally employ one-point statistical descrip- 
tion using Langevin type equation (a stochastic differen- 
tial equation of first order) for the dynamical variable, or 
the associated Fokker-Planck equation (a partial differen- 
tial equation) for one-point probability density function. 

With the choice of delta-correlated noises such mod- 
els fall into the class of Markovian models (no memory 
effects at small scales) allowing well established Fokker- 
Planck approximation. The consideration of finite-time 
correlated noises and the associated memory effects re- 
quires a deeper analysis which should be made separately 
in each particular case. The evolution equations are for- 
mulated and solved in the Lagrangian framework (the 
comoving frame), in a purely temporal treatment, with 
fluctuations being treated along the particle trajectory. 

Approximation of a short-time correlated noise by the 
zero-time correlated one is usually made due to the time 
scale hierarchy emerging from the general physical anal- 
ysis of the system and experimental data. Under the 
stationarity condition, a balance between the energy in- 
jected at large scales and the energy dissipated by viscous 
processes at small scales, one can try to solve the Fokker- 
Planck equation to find stationary probability density 
function of the acceleration, P(a). This function as well 
as the associated moments can then be compared with 
the experimental data on acceleration statistics. The 
Fokker-Planck approximation allows one to make a link 
between the dynamics and the statistical approach. In 
the case when stationary probability distribution can be 
found exactly one can make a further analysis without a 
dynamical reference, yet having a possibility to extract 
stationary time correlators. 

In contrast to the usual Brownian like motion, the 
fluid particle acceleration does not merely follow a ran- 
dom walk with a complete self-similarity at all scales. It 
was found to reveal a different, multiscale self-similarity, 
which can be seen from wide tails of a quasi-Gaussian 
distribution of the experimental probability distribution 
P(a). This requires a consideration of some specific 
Langevin type equations, which may include nonlinear 
terms, e.g., to account for turbulent viscosity effect, and 
an extension of the usual properties of model forces, ad- 
ditive and multiplicative noises. 

Specifically, the class of models represented by 
Eqs. Q-© is featured by consideration of the accel- 
eration evolution driven by the "forces" characterized 
by fluctuating drift coefficient 7 (or fluctuating intensity 
of multiplicative delta-correlated noise in a more gen- 
eral case) and/or fluctuating intensity er 2 of the additive 
noise. This was found to imply stationary distributions of 
the acceleration (or velocity increments in time, for finite 
time lags) of a quasi-Gaussian form with wide tails which 
are a classical signature of the turbulence intermittency, 
a phenomenon which developed turbulent flows exhibit 
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at small time scales. Earlier work on such type of models 
are due to Castaing, Gagne, and Hopfinger [lj], referred 
to as the Castaing model, in which a log-normal distribu- 
tion of fluctuating variance of intermittent variable was 
used without reference to a stochastic dynamical equa- 
tion. 

The difference from the wellknown class of stochastic 
models with delta-correlated Gaussian white multiplica- 
tive and additive noises which are also known to imply 
quasi-Gaussian stationary distributions with wide tails is 
that one supposes that intensities of the noises are not 
constant but fluctuate at a large time scale. We refer 
to the models with such Random Intensities of Noises as 
RIN models. 

This class of models introduces a two time-scale dy- 
namics, one associated to a delta-correlation of noises 
(modelling the smallest time scale under consideration, 
usually of the order of Kolmogorov time) and the other 
associated to variations of intensities of the noises, their 
possible coupling to each other, and other parameters 
assumed to occur at large time scales, up to a few La- 
grangian integral times. From a general point of view, 
one can assume a hierarchy of a number of characteristic 
time scales. However, in the present paper we simplify 
the consideration in order to make it more analytically 
tractable, in accord to the presence of two characteristic 
time scales in the Kolmogorov picture of fully developed 
turbulence. 

In the approximation of two time-scales, one can start 
with a Langevin type equation, derive the associated 
Fokker-Planck equation in Stratonovich or Ito formula- 
tions, and try to find a stationary solution of the Fokker- 
Planck equation, in which slowly fluctuating parameters 
are taken to be fixed. As the next step, one evaluates 
stochastic expectation of the resulting conditional prob- 
ability density function over the parameters with some 
distributions assigned to them. By this way one can ob- 
tain a stationary marginal probability density function 
as the main prediction of the model. 

The dynamical model J5J represents a particular sim- 
ple one-dimensional RIN model characterized by the 
presence of an additive noise (a short time scale) and 
fluctuating composite parameter j3 — j/cr 2 (a long time 
scale), where 7 is simply kinetic coefficient (a multiplica- 
tive noise is not present explicitly) and a 2 is the additive 
noise intensity. 

Two time-scale stochastic dynamics in describing the 
acceleration jointly with the velocity and position was 
used by Sawford [la ]. 

d t a = -(T i ; 1 +t- 1 )a + T [ ; 1 t- 1 u (10) 
+ ^2o-l{TZ l +tTi l )TZ% 1 L{t), 



where 
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d t x = u, 



_ 2a ^ 1 / 2 



(11) 



(12) 



are two time scales, Tl t v , Co, ao are Lagrangian 
structure constants, a 2 is the variance of the velocity 
distribution, and e is mean energy dissipation rate per 
unit mass. This model predicts Gaussian stationary dis- 
tributions for the acceleration and velocity reflecting un- 
corrected character of the fluctuations. An obvious ex- 
tension of this model is to replace e by stochastic energy 
dissipation rate e, and assume that it is log-normally dis- 
tributed in correspondence with the refined Kolmogorov 
1962 approach jl6j. 

Recently, an attempt to generalize the Sawford model 
to the case of fluctuating parameters has been made by 
A. Reynolds 0, with a very well agreement with the 
experimental results being achieved. 

The growing interest in studying Langevin type equa- 
tions to describe developed turbulence is motivated by 
the recent high precision Lagrangian experiments by 
Porta, Voth, Crawford, Alexander, and Bodenschatz 
[r?| . the new data by Crawford, Mordant, Bodenschatz, 
and Reynolds (the Taylor microscale Reynolds num- 
ber is R\ = 690, the normalized acceleration range 
is [—60,60] 9 a, the Kolmogorov timescale t v is re- 
solved) [lSf. Mordant, Delour, Leveque, Arneodo, and 
Pinton (R\ = 740, a e [-20,20], t v is not resolved) 
Mordant, Crawford, and Bodenschatz j^j, and direct 
numerical simulations of the Navier-Stokes equation by 
Kraichnan and Gotoh (R x = 380, a G [-150,150]) HJ; 
the classical Reynolds number is Re = R 2 /15. This gives 
an important information on the dynamics and new look 
to the intermittency in high-Rcynolds-numbcr fluid tur- 
bulence. 

Response characteristics of the polystyrene tracer par- 
ticle of about 46 /zm size and the precision in the ex- 
periments [rH ITsI ] allow to resolve about 1/20 of the 
Kolmogorov time and 1/20 of the Kolmogorov length in 
an R\ — 970 flow so that the acceleration can be re- 
ally resolved, and the particle follows rare violent events 
within 7% of the ideal value of acceleration even at 
the highest Reynolds number studied there. For lower 
Reynolds numbers the resolutions with respect to Kol- 
mogorov scales are even much higher. The collected 
statistics of about 1.7x10 s data points appeared to be 
sufficient to establish finitness of the fourth order mo- 
ment of the acceleration, (a 4 ). The acceleration values 
are obtained from the measured velocity increments in 
time by certain extrapolation to zero time increment, a 
procedure re quiri ng handling data points in the dissipa- 
tive scale fi7p2"o| . 

The stretched exponential fit with three parameters 
provides a very well agreement with the experimental 
data on the transverse acceleration a of the tracer particle 
in the Rx = 690 flow [HHI, 



P(a) = Cexp 



(l + \b ia /b 2 \ b3 )b 2 



(13) 



where b x = 0.513 ± 0.003, b 2 = 0.563 ± 0.02, and 
63 = 1.600 ± 0.003 are fit parameters, and C = 0.733 
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is a normalization constant. At large acceleration values 
the tails of the above P(a) decrease asymptotically as 
exp[— |a| - 4 ], that implies a convergence of the fourth or- 
der moment, (a 4 ) = a 4 P(a)da. The flatness factor 
of the distribution (|13|) which characterizes the widening 
of its tails (when compared with a Gaussian) is 



which should be compared with the flatness of the exper- 
imental curve, F = 55 ± 8 20]. We remind that for a 
Gaussian distribution F — 3. 

With constant /3, the Gaussian probability density 
function (|SJ| corresponds to the non-intermittent Kol- 
mogorov 1941 picture of fully developed turbulence, and 
agrees with the experimental statistics of components of 
velocity increments in time for large time scales, up to the 
integral time scale. However, it fails to describe observed 
Reynolds number dependent stretched exponential tails 
of the experimental acceleration probability density func- 
tion that correspond to anomalously high probabil- 
ities for the tracer particle to have extremely high accel- 
erations, bursts with dozens of root-mean-square (rms) 
acceleration, in the developed turbulent flow. Such a 
high probability of the extreme acceleration magnitudes 
is associated to the Lagrangian turbulence intermittency, 
which was found to be considerably stronger than the Eu- 
lerian one. Equivalently, one can say that it is related to 
an increase of the probability to have large velocity in- 
crements in time with a decrease of the time scale, down 
to the Kolmogorov time scale (a statistical viewpoint). 

In the Eulerian framework, the turbulence intermit- 
tency is usually understood differently, as an increase of 
the probability to have large longitudinal velocity differ- 
ences at short spatial scales, and studied through nonlin- 
earity in scaling exponents of velocity structure functions 
(a structural viewpoint). 

Intermittency of the stochastic energy dissipation rate 
is related to the dynamical intermittency of chaoticity in 
the system that makes a link between the Eulerian and 
Lagrangian intermittency to which we refer below. 

The averaging JQl of the Gaussian distribution © over 
randomly distributed f3, an evaluation of the stochastic 
expectation, was found to be a simple ad hoc procedure 
to obtain observable predictions, with one free parame- 
ter, that meet experimental statistical data on the accel- 
eration of the tracer particle. One can think of this as 
the averaging over a large time span for one tracer parti- 
cle, or as the averaging over an ensemble of tracer parti- 
cles, moving in the three-dimensional flow characterized 
by random spatially distributed domains with different 
values of (3. 

Physically one would like to know the processes un- 
derlying the random character of the model parameter 
(3. Due to the definition (0} the random character of (3 is 
attributed to a random character of the drift coefficient 
7 and/or the additive noise intensity a 2 . 



The distribution of (3 is not fixed uniquely by the the- 
ory so that a judicious choice of /(/?) makes a problem 
in the PJN model 0}-®. 

The primary aim of the present work is to provide a 
critical evaluation of the Langevin modeling approach 
by making a comparative analysis of different types of 
models on the basis of recent Lagrangian experimental 
data. 

The crucial point is to make a link between the 
Langevin type equation and the Navier-Stokes equation. 
This includes determination of statistical properties of 
stochastic terms and the functional form of deterministic 
terms, as well as their dependence on the parameters en- 
tering the Navier-Stokes equation justified for the inertial 
range of fully developed turbulence. Also, some extension 
of the stochastic equation may be required to account 
for dependence of the parameters on Lagrangian veloc- 
ity fluctuations, in the spirit of second-order stochastic 
models |15| and in correspondence to the Navier-Stokes 
equation as the pressure gradient term in the Eulerian 
framework can be expressed in terms of the velocity ow- 
ing to the incompessibility condition. Strong and nonlo- 
cal character of Lagrangian particle coupling due to pres- 
sure effects makes the main obstacle to derive turbulence 
statistics from the Navier-Stokes equation. 

The layout of the paper is as follows. 

In Sec. [n] we review some recent one-dimensional 
Langevin models of the developed turbulence. 

In Sec. Ill Al we outline implications of the Random In- 
tensity of Noises (RIN) models with the underlying chi- 
square fSec. IIIATt and log-normal (Sec. Ill A~2|l distribu- 
tions of/? HUGH We review results of a recent ap- 
proach jL| to specify f((3) which is based upon relating (3 
to velocity fluctuations u and using normal distribution 
of velocity fluctuations with zero mean f Sec. Ill A~3|) . This 
enables to reproduce chi-square and log-normal distribu- 
tions of (3 as particular cases. In general, this approach 
assumes that parameters of the model, such as the inten- 
sity of additive noise, depend on velocity fluctuations, 
in an agreement with the Heisenberg-Yaglom picture of 
developed turbulence. 

A nonlinear Langevin and the associated Fokker- 
Planck equations obtained by a direct requirement 
that the probability distribution satisfies some model- 
independent scaling relation have been recently proposed 
by Hnat, Chapman, and Rowlands [z| to describe the 
measured time series of the solar wind bulk plasma pa- 
rameters. We find this result relevant to fluid turbulence 
since it is based on a stochastic dynamical framework 
and leads to the stationary distribution with exponen- 
tially truncated power law tails, similar to that obtained 
in the above mentioned RIN models fSec. Ill B|l . 

The above one-dimensional Langevin toy models of La- 
grangian turbulence all suffer from the lack of physical in- 
terpretation, e.g., of short term dynamics, or small-scale 
and large-scale contributions, in the context of three- 
dimensional Navier-Stokes equation. 

The Navier-Stokes equation based approach to de- 
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scribe statistical properties of small scale velocity in- 
crements, both in the Eulerian and Lagrangian frames, 
was developed in much detail by Laval, Dubrulle, and 
Nazarenko [2j|; see also recent work [24|. This ap- 
proach introduces nonlocal interactions between well sep- 
arated large and small scales, elongated triads, and is 
referred to as the Rapid Distortion Theory (RDT) ap- 
proach. This approach is contrasted with the Gledzer- 
Ohkitani-Yamada shell model, in which interactions of a 
shell of wave numbers with only its nearest and next- 
nearest shells are taken into account. We outline re- 
sults of this approach and focus on the proposed one- 
dimensional Langevin model of Lagrangian turbulence 
to which we refer as Laval-Dubrulle-Nazarenko (LDN) 
model (Sec. Ill C|l . Particularly, we calculate exactly 
the probability density function of acceleration stemming 
from this model. 

In Sec. IIIII we make qualitative and quantitative com- 
parative analysis of the one-dimensional LDN and simple 
RIN models. 

In Sec. IIVI we study conditional probability density 
function P(a\u) taking the LDN model with delta- 
correlated noises and assuming that the additive noise 
intensity parameter a depends on the amplitude of ve- 
locity fluctuations u. 



II. ONE-DIMENSIONAL LANGEVIN MODELS 
OF THE LAGRANGIAN TURBULENCE 

In this Section, we outline results of some recent one- 
dimensional Langevin models of the Lagrangian fluid par- 
ticle acceleration in the developed turbulent flow. 



A. Simple RIN models 

1. The underlying chi-square distribution 

With the underlying gamma (chi-square) distribution 
of /? of order n (n — 1, 2, 3, . . .), 



f(P) 



r(f) V2/?o 



(3 2 1 exp 



n/3 
Wo 



(15) 



the resulting marginal probability density function 
with P(a\B) given by the Gaussian © is found in the 
form (cf. |6j) 



where 



P(a) = 



C 



c 



(a 2 +n//3 )(" +1 )/ 2: 
(n//3 )"/ 2 L(^) 



(16) 



(17) 



is normalization constant. With n — 3 (/?o = 3 for a unit 
variance) one obtains the normalized marginal distribu- 
tion in the following simple form, 



P(a) 



7r(a 2 + 1) 



2 ■ 



(18) 



a prediction of the chi-square model, with the Tsallis en- 
tropic index taken to be q = (n + 3)/(n + 1) = 3/2 due to 
the theoretical argument that the number of independent 
random variables at Kolmogorov scale is n = 3 for the 
three-dimensional flow . One can see that the resulting 
marginal distribution is characterized by power law tails 
that a priori lead to divergent higher moments. 

A Gaussian truncation of the power law tails naturally 
arises under the assumption that the parameter (3 con- 
tains a non- fluctuating part, which can be separated out 
as follows: (3/2 -> a~ 2 + 0/2 [HI This leads to the 
modified marginal distribution, 



P(a) 



Cexp[-a 2 /a 2 ] 
(a 2 +n/(3 Q )( n + i y 2: 



(19) 



where C is normalization constant and a c > is a free 
parameter which can be used for a fitting. Taking the 
theoretical value n = 3 and (3q = 3 as in the above case, 
one obtains that 



C 



2ai 



7r(a 2 - 2) exp[a^ 2 ](l - erf (ac 1 )) + 20ra c ' 



(20) 



where erf(x) denotes the error function. The distribu- 
tion fli^l at the fitted value a c = 39.0 is in a good agree- 
ment with the experimental probability density function 
P(a) 001. 

Note that at a c — ► oo (no constant part) the model 119|) 
covers the model (II 6|) . Within the framework of Tsallis 
nonextensive statistics, the parameter q — 1 measures a 
variance of fluctuations. For q — > 1 (no fluctuations), 
Eq. 1|19|) reduces to a Gaussian distribution, which meets 
the experimental data for temporal velocity increments 
at the integral time scale. 

A comparison of the chi-square model l|18|l and chi- 
square Gaussian model (|19f) with the experimental data 
is shown in Figs. Q an d El One can see that both the 
distributions follow the experimental P{a) to a good ac- 
curacy (at least up to 3 standard deviations), although 
the tails of the chi-square model distribution departure 
to the experimental curve at big |o|. A major difference 
is seen from the contribution to the fourth order moment 
a 4 P(a) shown in Fig. [21 The chi-square model yields a 
qualitatively unsatisfactory behavior indicating a diver- 
gency of the predicted fourth order moment. In contrast, 
the chi-square Gaussian model is in a good qualitative 
agreement with the data, reproducing them well at small 
and large acceleration values although quantitatively it 
deviates at intermediate acceleration values and gives the 
flatness value F ~ 46.1 for a c — 39.0, as compared to the 
flatness value (|14l) . 
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FIG. 1: Acceleration probability density function P(a). Dots: 
experimental data at R\ = 690 by Crawford, Mordant, and 
Bodenschatz 18]. Dashed line: stretched exponential fit 1131 1. 
6i = 0.513, b 2 = 0.563, b 3 = 1.600, C = 0.733. Dot-dashed 
line: Beck chi-square model (I18L q = 3/2. Solid line: Chi- 
square Gaussian model (1191 . a c = 39.0, C = 0.637. x = 
a/(a 2 ) 1 ^ 2 denotes normalized acceleration. 



2. The underlying log-normal distribution 



With the underlying log-normal distribution of /?, 



f(/3) 



27TS/3 



exp 



\2 



m ' 

2s 2 



(21) 



the resulting marginal probability density function l|T)l 
with P(a\0) given by the Gaussian © was recently pro- 
posed to be HI 



1 r°° 



(In 



P ^2 



2s 2 



e 2 



/3a" 



. (22) 



where the only free parameter, s, can be used for a 
fitting, or derived from theoretical arguments, s 2 = 3 
(to = exp[s 2 /2] for a unit variance). This distribution 
is shown in Fig. |3| and was found to be in a good agree- 
ment with the Lagrangian experimental data by Porta et 
al. Uffl, the new data by Crawford et al. Mordant et 
al. |l9j . and direct numerical simulations of the Navier- 
Stokes equation (DNS) by Kraichnan and Gotoh pl| . 



FIG. 2: Contribution to the fourth order moment a 4 P(a). 
Top panel: a linear plot, bottom panel: a log-log plot. Nota- 
tion is the same as in Fig. 



However, the central part of the distribution shown 
in the bottom panel of Fig. [3] reveals greater inaccuracy 
of the log-normal model (P(0) ~ 1.23) as compared with 
that of both the chi-square and chi-square Gaussian mod- 
els (P(0) ~ 0.65) which are almost not distinguishable in 
the region ^/(a 2 ) 1 / 2 < 4 (the bottom panel of Fig. QJ; 
see also recent work by Gotoh and Kraichnan 0] . This is 
the main failing of the log-normal model 1221) for s 2 = 3.0 
although the predicted distribution follows the measured 
low probability tails, which are related to the Lagrangian 
intermittency, to a good accuracy. The central region 
of the experimental curve H13|) (P(0) ~ 0.73) contains 
most weight of the experimental distribution and is the 
most accurate part of it, with the relative uncertainty of 
about 3% for \a\/(a 2 ) 1/2 < 10 and more than 40% for 
|a|/(a 2 ) 1/2 >40 IS. 

The distribution (|22|l is characterized by a bit big- 
ger flatness value, F = 3exp[s 2 ] ~ 60.3 for s 2 = 3, as 
compared to the flatness value l|14l) which is neverthe- 
less acceptable from the experimental point of view. The 
peaks of the contribution to fourth order moment shown 
in Fig.^ldo not match that of the experimental curve for 
the R\ — 690 flow. We note that the best fit is achieved 
for s 2 close to the theoretical value 3 but this does not 
significantly improve overlapping of the peaks with the 
data points. 

One naturally expects that a better correspondence to 
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FIG. 3: Acceleration probability density function P(a). Dots: 
experimental data at R\ = 690 by Crawford, Mordant, and 
Bodenschatz 18]. Dashed line: stretched exponential fit 1131 1. 
6i = 0.513, b 2 = 0.563, b 3 = 1.600, C = 0.733. Dot-dashed 
line: Beck log-normal model i'2'21 . s — 3.0. Solid line: Cas- 
taing log-normal model 123H . so = 0.625. x = a/la 2 ) 1 ' 2 de- 
notes normalized acceleration. 



the experiment may be achieved by an accounting for 
small scale interactions via turbulent viscosity (certain 
nonlinearity in the first term of the r.h.s. of Eq. P|l) as 
it implies a damping of the large events, i.e., less pro- 
nounced enhancement of the tails of P(a). 

It should be noted that the idea to describe turbulence 
intermittency via averaging of the Gaussian distribution 
over log-normally distributed variance of some intermit- 
tent variable was proposed long time ago by Castaing, 
Gagne, and Hopfinger (l4| . 



P(x) 



1 



2ttsq 



dee- 



' exp 



2sl 



-x 2 /(26 2 ) 



(23) 

where a; is a variable under study. Below, we apply this 
model to the Lagrangian acceleration, x = a. 

In technical terms, the difference from the Castaing 
log-normal model is that in Eq. (|22[1 the inverse square 
of the variance, (3 = 9~ 2 , is taken to be log-normally 
distributed. In essence, the models 1|22|) and (|23[) are 
of the same type, with different parameters assumed to 
be fluctuating at a large time scale and hence different 
resulting marginal distributions. 



FIG. 4: Contribution to the fourth order moment a 4 P(a). 
Top panel: a linear plot, bottom panel: a log-log plot. Nota- 
tion is the same as in Fig. [3] 



One can check that the change of variable, 9 = /3 -1 / 2 , 
in Eq. I|23ll leads to the density function different from 
that given by Eq. 



P(x) 



f 

Jo 



d/3(3 



-3/2 



exp 



(ln^) 2 

v rni ' 

2s\ 



(24) 



where we have denoted mi — and S\ = 2sq. There- 
fore, the distributions and are indeed not equiv- 
alent to each other, both being of a stretched exponential 
form. 

As to a comparison of the fits, we found that the fit of 
the Castaing log-normal model 123fl for the acceleration, 
with the fitted value s = 0.625 (m = exp [si/ 2} for unit 
variance), is of a considerably less quality as one can see 
from Fig. [3] and, more clearly, from Fig. 0] Positions 
of the peaks of a P(a) are approximately the same for 
both the models, namely, \a\j '(a 2 ) 1 / 2 ~ 8 as compared to 
^/(a 2 ) 1 / 2 ~ 10.2 for the experimental curve. 

We conclude this Section with the following remark. 
The Langevin model of the type Fokker-Planck ap- 
proximation of the type ©, and the underlying lognor- 
mal distribution l|21ll within the Castaing approach were 
recently used by Hnat, Chapman, and Rowlands j2^ to 
describe intermittency and scaling of the solar wind bulk 
plasma parameters. 
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3. The underlying Gaussian distribution of velocity 
fluctuations 

The problem of selecting appropriate distribution of 
the parameter /3 among possible ones was recently ad- 
dressed in Ref. 0. A specific model based on the as- 
sumption that the velocity fluctuation u follows normal 
distribution with zero mean and variance s was devel- 
oped. The result is that a class of underlying distribu- 
tions of P can be encoded in the function = P(u), and 
the marginal distribution is found to be 



P{a) = C{s) dj3 P(a\P) exp 
Jo 



HP)? 

2s 2 



du 
dp 



(25) 



where u(P) is the inverse function. Note that only an ab- 
solute value of u contributes this probability distribution. 
Particularly the exponential dependence, 



P(u) = exp[±it], 



(26) 



features the log-normal distribution of P so that Eq. 
leads to Eq. {221 used in Ref. 0, while the x 2 distribution 
of order one is recovered with 



P(u) 



(27) 



In general, this model is relevant when P(u) is a mono- 
tonic Borel function of the stochastic variable u mapping 
[—00,00] 9 u to [0,oo] 3 p, and allows one to rule out 
some ad hoc distributions of P as well as to make appro- 
priate generalizations of both x 2 and log-normal distri- 
butions of the parameter p. 

Below we develop a possible dynamical foundation of 
the above model. 

The stationary distribution (|25|1 with P(u) — exp[±w] 
can be associated to the Langevin equation of the 
form El, 



d t a = ~fF{a)+e u L{t), 



(28) 



where we denote lo — ^fu/2, u follows Gaussian distribu- 
tion with zero mean, and we take 7 = const to simplify 
the consideration. 

In general, we adopt a viewpoint that statistical prop- 
erties of the acceleration a are associated to velocity 
fluctuations statistics due to the wellknown Heisenberg- 
Yaglom theory. This theory predicts the following scaling 
of a component of the acceleration variance, 



= anM 9/2 */- 1/2 ir 3/2 , 



(29) 



where u is the rms velocity, L is the integral scale length, 
and ao is the Kolmogorov constant. This longstanding 
universal u 9 ^ 2 scaling was confirmed by the recent La- 
grangian experiments |17| to a very high accuracy, for 
about seven orders of magnitude in the acceleration vari- 
ance, or two orders of the rms velocity, at R\ > 500. 
At lower Reynolds numbers, R\ < 500, it appeared 
that the Heisenberg-Yaglom scaling significantly deviates 



from the experimental data due to the emerging depen- 
dence of ao on R\ (the acceleration is increasingly cou- 
pled to large scales of the flow at low Reynolds numbers). 

We note also that in the Sawford model the Langevin 
equation (|10|1 for a includes velocity fluctuations u and 
the variance of the velocity distribution. 

Below we outline a relationship of the model (|28|l to 
some recent approaches in studying the intermittency. 

(i) The form of the last term in Eq. (|28|l . in which ui 
can be viewed as a Gaussian process, uj — tu(t), indepen- 
dent on the white noise L(t), strikingly resembles that 
involved in the recently developed log-infinitely divisible 
multifractal random walks (MRW) model by Muzy and 
Bacry 25], a continuous extension of discrete cascades. 

(ii) The driving force amplitude of the form e u ^\ with 
an ultraslow decaying correlation function, (Lo(t)cu(t + 
r)} t — — Aq lnfr/Tj,], r < Tl, in Langevin type equa- 
tion has been recently considered by Mordant et al. [ljj ; 
Tj, stands for the Lagrangian integral time. The results 
of this model have been found in a very good agreement 
with the experimentally observed very slow decay of the 
equal-position time autocorrelation of the fluid particle 
velocity increment magnitude in time, |A r iti|, for each 
component (very much like to MRW model) attributed 
to the intermittency of Lagrangian trajectories in the de- 
veloped turbulent flow. Also, very slow decay was ob- 
served for the cross correlation of the magnitudes of the 
acceleration components. Both the dynamical correla- 
tions were found to vanish only for r > 3Tl, while the 
dynamical correlations of the full signed entities, A T Ui, 
decay rapidly, the autocorrelation functions cross zero at 
about r ~ 0.06T7 and the cross correlation functions are 
approximately zero. 

We note that the fitted value of the above intermit- 
tency parameter Ao (A§ = 0.115 ± 0.01) is very close to 
1/s 2 = 1/3, with the value s 2 = 3 in Eq. |22J) interpreted 
as the number of independent random variables in three- 
dimensional space at the Kolmogorov scale. If this is 
not due to a coincidence, the intermittency parameter 
Ao approaches simply the inverse of the effective space 
dimension number d, 



Ao = l/d, 



(30) 



d = 3, for high-Reynolds-numbcr turbulent flows. 

The above connection to the MRW model and very 
slow decay of the correlations of the absolute values of 
acceleration components indicate relevance of the spe- 
cific representation l|28|l . with very slow varying \u\, in 
the description of intermittency. In fact, due to the ex- 
periments 26] the Lagrangian velocity autocorrelatation 
function (u(t)u(t + T))t/(u 2 ) decays almost exponentially 
but very slowly, to vanish only for r > 317,, where the 
integral time scale Tl = 2.2 x 10~ 2 s, which is two orders 
of magnitude bigger than the Kolmogorov time scale, 
t„ = 2.0 x 10~ 4 s; i? A = 740 and the mean velocity 
is about 10% of the rms velocity. 

(iii) Due to the wellknown Kolmogorov power law scal- 
ing relationship between e and u, the representation l|28|) 
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can be thought of as the result of use of the relation 
In/? ~ me, with lne being normally distributed due to 
the refined Kolmogorov 1962 theory. Here, e denotes the 
stochastic energy dissipation rate per unit mass treated 
in the Lagrangian framework. From this point of view, 
one can identify 

cu = g In e, (31) 

where g is a constant. This means that the stochastic 
dynamics of the logarithm of the energy dissipation is 
independent, and it influences the acceleration dynam- 
ics specifically through the intensity of driving stochastic 
force in Eq. 1)28(1 . Stationary normal distribution of lo 
can be in turn derived from the Fokker-Planck equation 
associated to the Langevin equation of a linear form, 

d t uj = g +g\U) + g 2 L(t), (32) 

where gi are constants. This equation is in an agreement 
with the recent results of Eulerian (hotwire anemometer) 
study of the interaction between velocity increments and 
normalized e nerg y dissipation rate by Renner, Peinke, 
and Friedrich [23 • Particularly, they found that an expo- 
nential dependence of the diffusion coefficient on the log- 
arithmic energy dissipation in the Fokker-Planck equa- 
tion for the velocity increments in space is in a very good 
agreement with the experimental data. We note that 
Eq. (|32|) does not imply a logarithmic decay of the La- 
grangian correlation function (tu(t)uj(t + r)) proposed in 
Ref. [13 • This may be attributed to the wellknown differ- 
ence between the Eulerian (fixed probe) and Lagrangian 
(trajectory) frameworks. 

(iv) For the choice (3(u) = exp[it] corresponding to the 
log-normal distribution of /?, using Eq. (|28|l one can de- 
rive the stationary probability density function of the 
form [l2| 

/oo 
du C{u) exppn[ff(u)] - eV/2], (33) 
-00 

where g(u) is a probability density function of u. Hence 
the joint probability density function can be written 

P(a, u) = C{u) exp[ln[ ff (u)] - e u a 2 /2] . (34) 

Such a form of the distribution, containing specifically 
the double exponent, resembles the "universal" distribu- 
tion of fluctuations (Gumbel function), 

P(x) = c exp[ci(y - e y )}, y = c 2 (x~c 3 ), (35) 

where Ci are constant, recently considered by Chapman, 
Rowlands, and Watkins |2g (see also references therein) 
following the work by Portelli, Holdsworth, and Pin- 
ton (2^. They used an apparently different approach 
(not related to a Langevin type equation for the accel- 
eration studied in our paper) based on the multifractal 
type energy cascade and \ 2 or log-normal (K62 theory) 
underlying distribution for the energy dissipation rate at 



fixed level. They pointed out a good agreement of such 
P(x) with experimental data, where x denotes a fluc- 
tuating entity observed in a variety of model correlated 
systems, such as turbulence, forest fires, and sandpiles. 
The result of this approach meets ours and we consider 
it as an alternative way to derive the characteristic prob- 
ability measure of fluctuations; with g(u) taken to be a 
X 2 (respectively, Gaussian) density function, one obtains, 
up to a pre-exponential factor and constants, P(u) ~ 
exp[— u~ exp[it]] (respectively, P(u) ~ exp[— u 2 — exp[u]]). 
Thus we conclude that the "universal" distribution 135|) 
can be derived also within the general framework pro- 
posed in Ref. [12J that reflects a " universal" character of 
the underlying \ 2 distribution Q. 

We note that although successful in describing the ob- 
served statistics of Lagrangian acceleration, with a few 
simple hypotheses and one fitting parameter, the one- 
dimensional Langevin RIN models GJ-© and ffi fy -pg jl 
suffer from the lack of physical interpretation in the con- 
text of three-dimensional Navier-Stokes equation. 

In summary, we have presented a class of models 
(|28|l using the basic assumption that the parameter j3 de- 
pends on normally distributed velocity fluctuations. This 
class has been found to incorporate the previous RIN 
models in a unified way, with the dependence (3(u) re- 
quired to be a (monotonic) Borel function of the stochas- 
tic variable u. 



B. Hnat-Chapman-Rowlands model 

Another interesting model developed recently by Hnat, 
Chapman, and Rowlands j2^ to describe the observed 
time series of the solar wind bulk plasma parameters is 
based on the construction of Fokker-Planck equation for 
which the probability density function obeys the follow- 
ing one-parametric model-independent rescaling, 

P(x, t) = t- aa P s (xt- ao ), (36) 

where x denotes fluctuating plasma parameter and a 
is the scaling index. The value ao = 1/2 corresponds 
to a self-similar Brownian walk with Gaussian probabil- 
ity density functions at all time scales. The fitted value 
ao = 0.41 corresponds to a single non-Gaussian distri- 
bution P s (x s ), to which the observed distributions of 
some four plasma parameters collapse under the scaling, 

The Langevin equation of this model assumes only ad- 
ditive noise, and in such an anzatz it was found to be 

d t x = Dx{x)+D 2 {x)ri{t), (37) 

where D\ and D 2 are of the form 



Di{x) = 




D 2 (x) = (6 (1 - X ) - ao)* 1- " " 1 , (39) 
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do, bo are constants, and 2Dq is intensity of the delta- 
correlated Gaussian- white additive noise rj(t), (r)(t)) — 
0. By construction, this specified form of the dynamical 
equation ensures that the corresponding Fokker-Planck 
equation, 

d t P{x,t) =d x [a x 1 - a o 1 p + b x 2 - a o X d x P], (40) 

has the general solution P(x, t), which exhibits the scal- 
ing The fitted values are ao/&o = 2, bo = 10, 
and a^ 1 = 2.44. The rescaled distribution P s (x s ) corre- 
sponding to Eq. (|40Jl is characterized by power law tails 
truncated by stretched exponential with a good fit to the 
tails of the experimental distribution, but it diverges at 
the origin, x s — ► 0. 

To sum up, we point out that this diffusion model uses 
the generalized self-similarity principle resembling that 
used in the Eulerian description of the energy cascade 
in developed three-dimensional fluid turbulence, and ap- 
pears to be valid only asymptotically for large values of 
the variable, with the fitted parameter value being about 
a = 0.41. 

C. Laval-Dubrulle-Nazarenko model 

The Navier-Stokes equation based approach to de- 
scribe statistical properties of small scale velocity incre- 
ments, both in the Eulerian and Lagrangian framework, 
was develo ped in much detail by Laval, Dubrulle, and 
Nazarenko |23j | ; see also recent work |24| . This approach 
is based on featuring nonlocal interactions between well 
separated large and small scales, elongated triads, and 
referred to as the Rapid Distortion Theory (RDT) ap- 
proach. Decomposition of velocities into large and small 
scale parts was made by introducing a certain spatial 
filter of a cut off type. Within the framework of this 
approach, a three-dimensional Langevin model of the de- 
veloped turbulence was proposed. 

In its one-dimensional version, the toy model of the 
Lagrangian turbulence naturally implies a nonlinear 
Langevin type equation for a component of the small 
scale velocity increments in time (in the zero time-scale 
limit it corresponds to the acceleration a of fluid parti- 
cle) H, 

d t a= {£-is t k 2 )a + <j ± . (41) 

This equation is a Lagrangian description in the scale 
space, in the reference frame comoving with a wave num- 
ber packet. Here, 

v t = yJv%+B*a?/k 2 (42) 

stands for the turbulent viscosity introduced to describe 
small scale interactions, Vq is kinematic viscosity, B is 
constant, k is wave number (d t k — — k(0) = k , 
to model the RDT stretching effect in one-dimensional 
case), £ and a± are multiplicative and additive noises 



associated to the velocity derivative tensor and forcing 
of small scales by large scales (the energy transfer from 
large to small scales), respectively. 

We refer to the model (|41|) as the one-dimensional 
Laval-Dubrullc-Nazarenko (LDN) model of the devel- 
oped Lagrangian turbulence. This toy model can also 
be viewed as a passive scalar in a compressible one- 
dimensional flow. 

The noises, one-dimensional versions of which appear 
in Eq. (|4T)l . are projections related to the large-scale ve- 
locities, IL, and small-scale velocities, m, of the flow as 
follows 

Ik -> ^ 

£ = V (-(fcE/)-E/), (43) 
k -> 

& ± =&-—(ka), (44) 



(Ti = dj {Ui Uj - Ui Uj + Uj Ui - U jUi ) , (45) 

where the hat denotes Gabor transformation [30j, and 
the bar stands for the spatial cut off retaining the large 
scale part. One can see that the noises are related to 
the velocity fluctuations and the additive noise contains 
interaction terms between the large and small scale dy- 
namics. 

This gives a support to the idea that the intcrmittency 
is caused also by some nonlocal interactions within the 
inertial range and not merely by small scales. We remark 
that one would also like to know the role of the dissipative 
scale in this integrated picture. 

Noisy character of the entities Q43[) and l|44|l may not 
be seen as a consequence of the Navier-Stokes equation, 
which does not contain external random forces at the 
characteristic time scale. In the RDT approach, £ and 
a±_ are treated as independent stochastic processes en- 
tering the small scale dynamics (|41ll owing to the fact 
that the large-scale dynamics is weakly affected by small 
scales (that corresponds to a direct energy cascade in the 
three-dimensional flow) and thus can be viewed as a given 
noise. 

The relations l|43|l and l|44|l can be used to trace back 
the origin of the multiplicative and additive noises enter- 
ing various surrogate Langevin models of the developed 
turbulence, and to provide important information on the 
dynamics underlying the intermittency. 

Statistical properties of all the components of the 
noises £ and a± were studied numerically for decaying 
turbulence and reveal rich and complex behavior, in the 
laboratory frame. 

As a first step in the one-dimensional case, these noises 
were modelled in the Lagrangian frame by the coupled 
delta-correlated Gaussian white noises |2^| , 

(£(*)> = 0, mm) = 2D5(t-n 

(a ± (t)) = 0, {<T ± (t)a ± (t')) = 2a6(t - i'), (46) 
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where D, a, and A are parameters depending on scale via 
k . Stationary solution of the Fokker-Planck equation 
associated to Eq. (ITTJl with the noises iffify . 

d t P{a,t) = d a {v t k 2 P) + Dd a {ad a aP) - \d a {ad a P) 

-Xd 2 a (aP)+ad 2 a P, (47) 



is given by [2S 



P{a) = Cexp 



dy- 



-v t k 2 y - Dy + X 



Dy 2 - 2Xy + a 



(48) 



where C is a normalization constant and six parameters 
can be used to fit the experimental data. This model 
specifies the one-dimensional LDN model (|41[) , and we re- 
fer to this model as the LDN model with delta-correlated 
noises (dLDN model). 

The Langevin equation containing both the delta- 
correlated Gaussian white multiplicative and additive 
noises was studied in detail by Nakao [3l| (see also ref- 
erences therein) by using the associated Fokker-Planck 
equation. The dLDN model l|41ll extends Nakao's set up 
by incorporating two new features: (i) the nonlinearity 
controlled by B in Eq. I|42l) and (ii) the coupling of the 
noises controlled by A in Eq. I|46|l . 

It is interesting to note that the RDT approach qual- 
itativel y re sembles the model studied by Kuramoto and 
Nakao j32j, a system of spatially distributed chaotic el- 
ements driven by a field produced by nonlocal coupling, 
which is spatially long-wave and temporally irregular. 
Such systems, in which the multiplicative noise is the lo- 
cal Lyapunov exponent fluctuating randomly due to the 
chaotic motion of the elements, show power-law correla- 
tions, intermittency, and structure functions similar to 
that of the developed fluid turbulence. 

In the next Section, we make a comparison of the RIN 
model CD - © with the Laval-Dubrulle-Nazarenko model 
(|41|l . as well as its particular case, the dLDN model (|48ll . 



III. COMPARISON OF THE SIMPLE RIN AND 
LDN MODELS 

A. A qualitative comparison 

A direct comparison of the Langevin equations @ and 
(|41|l of the two models suggests the following evident 
identifications: 



"fF(a) — (£ - v t k 2 )a, oL = a±. 



(49) 



Hence the additive noises can be made identical to each 
other by putting a 2 = a. Further, in the case of a linear 
drift force, F(a) = —a, and constant viscosity, v% = vq, 
we can identify the remaining parameters, 7 = v$k 2 — £, 
so that we get 



/3 = 1 /a 2 = (v k 2 - 0/a. 



(50) 



This relation implies that the parameter (3 can be viewed 
as a stochastic variable with a nonzero mean due to the 
stochastic nature of £ assumed in the LDN model. This 
is in agreement with the simple RIN model, the defining 
feature of which is just that the fluctuating part of (3 
follows some statistical distribution. 

In the dLDN model |gT}-J2H||, both the additive and 
multiplicative noises are taken delta- correlated due to 
Eq. (|46|) . This is in a sharp contrast to the assumption 
that P can be taken constant to derive the stationary so- 
lution (JSJ) which is in the foundation of the simple RIN 
model. More precisely, the solution in the form © can 
be obtained as the lowest order approximation if (3 is 
slow varying in time as compared to a typical time scale 
associated to the additive noise L(t) (the adiabatic ap- 
proximation) . This suggests that the multiplicative noise 
£ should be taken sufficiently slow varying stochastic vari- 
able, to meet the anzatz used in the RIN model. 

The detailed numerical analysis of the noises (2j| for 
the turbulent flow at relatively low Reynolds numbers, 
57 < R\ < 80, shows that the autocorrelation of the 
multiplicative noise £ decays much slower (by about one 
order of magnitude) than that of the additive noise a± . 
Hence a typical time scale at which £ varies, r^, is con- 
siderably bigger than that, t ct , of a±. Also, the cross- 
correlation between the two noises was found to be rather 
weak, A <C D and A <C a, by about two orders of magni- 
tude in the longitudinal case, and A = in the transverse 
case. Altogether this allows one to introduce the time 
scale hierarchy, ^ T a , and to decouple the noises, i.e. 
to put A = 0, that justifies the adiabatic approximation 
and the one-dimensional RIN model. 

The presence of the longtime correlated amplitude, 
e u ^\ and the short-time correlated directional part, L(t), 
of the stochastic driving force in Langevin type equation 
considered by Mordant et al. jljj also supports the above 
adiabatic approximation (two well separated time scales 
in the single additive stochastic force, in the Lagrangian 
framework). As usual, the delta-correlated noise origi- 
nates from taking the limit of zero correlation time in a 
system with the smallest finite noise correlation time. 

On the contrary, in the dLDN model one assumes 
the approximation of comparable time scales, ~ T a , 
and retains the coupling parameter A relating small scale 
stretching and vorticity and responsible for the skewness, 
which is however quite small in homogeneous isotropic 
turbulent flows. 

The use of the constant turbulent viscosity, v t = i/q, 
makes a good approximation in describing intermittency 
corrections since both the constant and turbulent viscos- 
ity were found to produce corrections which are of the 
same level as the DNS result [23j. In the physical con- 
text, this means that the small scale interactions are not 
of much important in the dynamics underlying the in- 
termittency. This justifies the use of the approximation 
of linear forcing, F(a) = —a, in the simple RIN model. 
We note that this is also in an agreement with both the 
experimental results for the Lagrangian velocity auto- 
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correlation function by Mordant, Metz, Michel, and Pin- 
ton |26| , and the recent experimental Eulerian results for 
the spatial velocity increments by Renner, Peinke, and 
Friedrich [27| . 

Alternatively, one can consider a more general RIN 
model characterized by the presence of delta-correlated 
Gaussian white additive and multiplicative noises and 
fluctuating intensities of both the noises. This will lead 
to the model similar to the dLDN model (|41|l in which 
the noise intensities D and a, and the coupling parameter 
A are assumed to fluctuate at a large time scale. 

In summary, we found that the one-dimensional RIN 
model 0-© can be viewed as a particular case of the 
one-dimensional LDN model (|41|) of turbulence which is 
based on the RDT approach by Laval, Dubrulle, and 
Nazarenko t 23] . It should be stressed that while both the 
toy models assume introduction of some external statis- 
tics — the correlator of L{t) and the distribution f({3) 
in Eq. and the correlators of £ and a± in Eq. (|41|l — 
the LDN model is characterized by a solid foundation and 
reveals a rich structure as compared to the RIN model. 

In the first approximation, i.e., A = 0, ft = is®, and 
T£ ^> T a , the class of RIN models is in a quite good qual- 
itative correspondence with the LDN model l|41|) and dif- 
fers from the specific dLDN model (|4"T)l - l|4"g|) by the only 
fact that in the latter one assumes ~ T a and intro- 
duces a delta-correlated multiplicative noise. Hence the 
different resulting probability density functions for the 
acceleration of fluid particle in the developed turbulent 
flow, Eqs. (Hnj-122) and |@BJl, respectively. 



B. A quantitative comparison 

With the above result of the qualitative comparison, we 
are led to make a more detailed, quantitative comparison 
of the dLDN model gTJ-(EEJ> and the simple RIN model 
Q-© with the underlying \ 2 or log-normal distribu- 
tion of /3, in order to determine which approximation, 

~ T a or ^> r a , is better when used to describe the 
Lagrangian statistical properties of the developed turbu- 
lent flow. We take the recent high precision Lagrangian 
experimental data 0, on statistics of fluid particle 
acceleration in the developed turbulent flow as a testbed. 
Actually we follow the remark made in Ref. [2^ that the 
delta approximation of £ is debatable and the perfor- 
mance of such a model should be further examined in 
the future. 

In Ref. [2i|, explicit analytic evaluation of the distri- 
bution (|48|l is given for the particular case v t = i>o, while 
the general case is treated in terms of d\nP(a)/da when 
fitting to the numerical RDT data (—4 < a < 4). In or- 
der to make fits to the experimental probability density 
function P(a) and to the contribution to the fourth order 
moment, a 4 P(a), covering wide range of the normalized 
acceleration, —60 < a < 60, one needs in an analytic or 
numerical evaluation of the r.h.s. of Eq. (|48|l . To this 
end, we have calculated exactly the integral appearing 



in the dLDN probability density function 1481) (see Ap- 
pendix A). 

The x 2 and log-normal distribution-based probability 
density functions (|19fl and (|22|l are both realizations of 
the RIN model, and contain one fitting parameter, a c and 
s, respectively. The result of comparison of fitting quali- 
ties of these functions ^l| 1 with a c = 39.0 and s = 3.0, is 
that the probability density function l|22|) provides a bet- 
ter fit to the experimental data 0| on low-probability 
tails and the contribution to the kurtosis summarizing 
the peakedness of distribution. However, since the in- 
tegral in Eq. (1221) can not be evaluated analytically we 
will use the distribution (|19fl . which provides a better fit 
to the central region, when dealing with analytic expres- 
sions. 

The dLDN probability density function (|48|l contains 
six parameters which can be used for a fitting, the mul- 
tiplicative noise intensity D, the additive noise intensity 
a, the coupling A between the multiplicative and additive 
noises, the turbulent viscosity parameter B, the param- 
eter fo, and the wave number parameter k. 

The parameter k. For a fitting, we can put k = 1 
without loss of generality since it can be absorbed by the 
redefinition of the parameters Vq and B, 

v Q k 2 v , Bk B. (51) 

The parameter a. The structure of the r.h.s. of 
Eq. H48fl is such that only four parameters out of five can 
be used for a fitting. For example, one can put a = 1 
without loss of generality by using the following redefini- 
tions, 

v n ja -> 2/ , B/a B, D/a D, A/a A. (52) 

Alternatively, one can put D = 1 provided the redefini- 
tions 

v /D v , B/D B, a/D -> a, X/D A. (53) 

The parameter A. Due to Eq. i|A5(l we have c = 
— i\J Da. — A 2 , which is pure imaginary for Da > A 2 and 
real for Da < A 2 . For c = 0, i.e. Da = A 2 , the integral 
(IA3I) is finite since divergent F(c) and F(—c) defined by 
Eq. i|A4() cancel each other. Since the parameter A mea- 
suring the coupling between the noises is assumed to be 
much smaller than both the noise intensities D and a [2^| , 
we put Da > A 2 in our subsequent analysis. Moreover, 
the parameter A responsible for the skewness can be set 
to zero since we will be interested, as a first step, in sta- 
tistically isotropic and homogeneous turbulent flows, for 
which the experimental distribution P(a) exhibits very 
small skewness }l7j . 

Thus, we can use three redefined free parameters, iso, 
B, and D, for a fitting, with k = 1, a = 1, and A = 0. 
However, we shall keep k and a in an explicit way in the 
formulas below, to provide a general representation. 

We start by considering two important particular cases 
of the dLDN probability density function (|48|) : a constant 
viscosity, v t = ^0, and a dominating turbulent viscosity, 
v t = B\a\/k. 
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1. Constant viscosity 

At A = (symmetric case) and B = 0, i.e. con- 
stant viscosity v t = v Q , using Eq. (|A1|) in Eq. ((4*S|) we 
get (cf. El) 



P{a) = C(Da 2 + a )-(i+"ofc 2 /o)/2 7 



(54) 



where C is normalization constant. This distribution is 
of a power law type and we can compare it with the re- 
sult I jl 9(1 , which contains a Gaussian truncation of similar 
power law tails. 

We note that with the identifications, 

D/a = 2(q-l), (1 + u Q k 2 /D)/2 = l/(q - 1), (55) 

the distribution 154(1 coincides with that obtained in the 
context of generalized statistics with the underlying x 2 
distribution 6]. Particularly, for q = 3/2 (n = 3, j3o = 3) 
used there, it follows that D/a = 1 and v k 2 = 3. 

It is highly remarkable to note that the two different 
approaches yield stationary distribution of exactly the 
same power law form, for certain identification of the pa- 
rameters. Namely, the Gaussian white delta-correlated 
multiplicative and additive noises with constant intensi- 
ties and a linear drift term imply P(a) of the same form 
as that obtained in the RIN model with \ 2 distributed 
/?, the ratio of the drift coefficient to the intensity of the 
Gaussian white delta-correlated additive noise. It follows 
that the effect of x 2 distributed mimics the presence of 
the multiplicative noise, and vice versa, in this particular 
case. 

The power law distribution ((54(1 can be used to get 
a good fit of the Lagrangian experimental P(a) data for 
small accelerations, e.g., with the normalized values rang- 
ing from —10 to 10, but in contrast to the Gaussian trun- 
cated one 1(19(1 it exhibits strong deviations for large a, 
and for (1 + v k 2 /D)/2 < 2 leads to a divergent fourth 
order moment, which is known to be finite [Tol. ITU . 

Introducing the noise intensity ratio parameter, 



and denoting 



« = -(! + v k 2 /D)/2 



(56) 



(57) 



we can rewrite the normalized distribution ((54|) as follows 

(cf. 0) 



P(a) = 



(l + b 2 a 2 Y 



2 2 F 1 (- K ;±;§;-&2)' 



(58) 



where %Fi is the hypergeometric function. In accord to 
the analysis made by Nakao |3lj | , for small additive noise 
intensity, i.e., at b 3> 1, this distribution exhibits a pro- 
nounced plateau near the origin, and the nth order mo- 
ments, truncatedby reflective walls at some fixed \a\, be- 
have as a power of b, 

(a n ) ~ b- Uok2/D , for n > v k 2 /D, (59) 
(a") ~ ft-", for n < v k 2 /D, (60) 



where n > 0. Thus, the truncated moments behave as 
(a n ) ~G + Gib-* (n) , (61) 

where Go^ are some constants and the function H(n) is 
zero at v^k 2 = and monotonically saturates to n at big 
vq k 2 . It should be stressed that such a behavior of the 
moments for small additive noise intensity is not specific 
to the distribution ((58(1 since it gives divergent moments 
but arises after some truncation of it, for example, by 
means of reflective walls or nonlinearity. Particularly, a 
truncation of the power law tails of the distribution nat- 
urally arises when accounting for the turbulent viscosity 
to which we turn below. 

2. Dominating turbulent viscosity 

At A = (symmetric case), for the case of dominating 
turbulent viscosity, v t = B\a\/k, using Eq. ((A2() we get 
for positive and negative a, respectively, 

(-i e fBka/D±Bka 1/2 D-' i/2 &Ycta,n[(D/a) 1/2 a\ 

P ^ = (Da 2 + ayf 2 ' (62) 

where C is normalization constant. One can see that, 
as expected, the power law dependence is of a similar 
form as in Eq. 1(54(1 but it is exponentially truncated at 
big |a| owing to the turbulent viscosity term. This dis- 
tribution is similar to the Gaussian truncated one l(19|) 
but the truncation is of an exponential type and there 
is some symmetric enhancement of the tails supplied by 
the arctan term. 

Now we turn to the general case, which provides a 
link between the two particular cases, v t — i/q and v t = 
B\a\/k, considered above. 



3. The general symmetric case 

At A = (symmetric case), from Eqs. I|A5(I - (|A7(I we 
have 



c = — id,2, c\ = id\d,2, c,2 = kd\, 
where we have denoted 



(63) 



(64) 



d l = JD(Dk 2 v 2 - B 2 a), d 2 



Note that c is pure imaginary and the r.h.s. of l(A3ll is 
much simplified yielding a symmetric distribution with 
respect to a — > —a. The entity a defined by Eq. ((A7() 
may be either real (for Dk 2 v^ > B 2 a) or pure imaginary 
(for Dk 2 v$ < B 2 a). In particular, for the case of con- 
stant viscosity, v$ 3> B 2 , it is real while for the case of 
dominating turbulent part of the viscosity, B 2 3> f 2 , it 
is pure imaginary, provided that the intensities of noises, 
D and a, are of the same order of magnitude. These two 
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particular cases lead to different final expressions for the 
distribution, (|54|) and (|62[) . respectively, obtained above. 

In the general case, using Eqs. (|63|l and l|64JI in 
Eq. ijSHIl after some algebra we obtain the following ex- 
pression for the dLDN probability density function l|48|) . 
at A = 0, 



P(a) 



4D 5 (B A Daa 2 



Ce-^ k2 ' D 
{Da 2 + a) 1 / 2 

k z {Dkv^ + d^tY) 1 



dfk 2 (Da 2 + a) 



(65) 



where C is normalization constant, d\ is given in Eq. (16411 . 
and vt is given by Eq. 142(1 . 

It can be easily checked that Eq. (|65[1 reduces to 
Eq. I|54|) at B = 0, while to verify that it reduces to 
Eq. (|62[1 at v t = B\a\/k requires the use of the fact that 
d\ becomes pure imaginary, returning back to the log- 
arithmic representation due to Eq. (|A4|) . and the iden- 
tity ijXSJl. 

The distribution (|65fl is characterized by the power law 
tails, which are (i) exponentially truncated and (ii) en- 
hanced by the power law part of the numerator, with 
both the effects being solely related to the nonzero tur- 
bulent viscosity coefficient B responsible for a nonlinear 
small scale dynamics. 

We conclude that to provide an acceptable fit of the 
dLDN model prediction to the Lagrangian experimental 
data [l8l | nonlinear small scale interactions encoded in 
the turbulent viscosity v t are essential. 

Sample fit of the dLDN probability density function 
P(a) given by Eq. (|65|1 and contribution to the fourth or- 
der moment, a 4 P(a), are shown in Figs. |S] and H3 respec- 
tively. In the numerical fit, we have put, in accord to the 
redefinitions (|52|) , the wave number parameter k = 1 and 
the additive noise intensity parameter a — 1 in Eq. I|65|) 
and fitted the remaining three parameters, vq, D, and 
B. One can observe a good agreement with the exper- 
imental data. Particularly, the dLDN contribution to 
the kurtosis a A P{a) plotted in Fig. does peak at the 
same points as the experimental curve (positions of the 
peaks depend mainly on D). The central part of the 
dLDN distribution shown in the bottom panel of Fig. [5] 
fits the experiment to a higher accuracy as compared 
with the log-normal model i|22[l but yet departures from 
that of the experimental curve. This departure can be 
attributed to the approximation of delta-correlated mul- 
tiplicative noise used in the dLDN model (see discussion 
in Sec. IIIII above). 

Having the general form of the dLDN distribution eval- 
uated explicitly, Eq. (|65|) , one can derive higher accelera- 
tion moments, (a™), n = 2, 4, . . . The associated integrals 
are not analytically tractable and can be evaluated nu- 
merically. We will consider these in Sec. II VI below. 

In the most general case (A ^ 0) the resulting P(a) 
is given due to an exponential of the exact integral (|A3|) 
which we do not represent here for brevity. 
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FIG. 5: Acceleration probability density function P(a). Dots: 
experimental data at R\ = 690 by Crawford, Mordant, and 
Bodenschatz [lfijl . Dashed line: stretched exponential fit 1131 1. 
6i = 0.513, b 2 = 0.563, b 3 = 1.600, C = 0.733. Dot-dashed 
line: Beck log-normal model (1221 . s = 3.0. Solid line: Laval- 
Dubrulle-Nazarenko model (1651 . k = 1, a = 1, D = 1.130, 
B = 0.163, v = 2.631, C = 1.805. x = a/(a 2 ) 1/2 denotes 
normalized acceleration. 
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FIG. 6: Contribution to fourth order moment a 4 P(a). No- 
tation is the same as in Fig. 



To sum up, we have made an important step forward 
with the dLDN model by having calculated P(a) exactly. 
We have shown that the dLDN model is capable to re- 
produce the recent Lagrangian experimental data on the 
acceleration statistics to a good accuracy. Particularly, 
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we found that the predicted fourth order moment den- 
sity function does peak at the same value of acceleration, 
\a\/(a 2 ) x / 2 ~ 10.2, as the experimental curve, in con- 
trast to the predictions of the other considered stochastic 
models. The presence of the delta-correlated multiplica- 
tive noise and the nonlinearity (turbulent viscosity) in the 
model Langevin equation was found to be of much impor- 
tance. The considered RIN models provide less but yet 
acceptable accuracy of the low-probability tails despite 
they employ only one free parameter, which can be fixed 
by certain phenomenological arguments, as compared to 
the dLDN model, which contains four free parameters. 
However, we stress that in contrast to the LDN model 
the considered RIN models have a meager support from 
the turbulence dynamics. 



IV. CONDITIONAL PROBABILITY DENSITY 
FUNCTION OF THE ACCELERATION 

In the very recent paper [2(| a new experimental data 
on the conditional probability density function of the 
transverse acceleration, P(a\u), have been reported. In 
general, this representation is in agreement with the pro- 
posed idea that velocity fluctuations u are directly in- 
volved in the stochastic acceleration dynamics 0] rep- 
resented in Sec. Ill A 31 

Also, the observed conditional acceleration component 
variance has been found to be in a good agreement with 
the Sawford et al. scaling relation (see [2[j and references 
therein) , 

(a»~w 6 , (66) 

obtained to a leading order in the same component u (to 
be not confused with the rms velocity u = (w 2 ) 1 / 2 ). 

The experimental data reveal highly non-Gaussian, 
stretched exponential character of P(a\u), very similar 
to that of P(a), for fixed u ranging from zero up to three 
rms velocity U :2jj as opposed to the theoretical result 
that P(a\u) is a Gaussian in a due to the simple RIN 
model ©, with arbitrary (3 = (3(u), or due to the more 
general RIN model (1331) . Similarity between the exper- 
imental P(a\u) and P(a) suggests that they share the 
process underlying the fluctuations. 

Below we address this important problem within the 
framework of the RIN approach. 

The idea is that the stretched exponential form of the 
tails of the observed conditional distribution P(a\u) could 
be assigned solely to small time scales, while the marginal 
probability distribution P(a) is developed from P(a\u) at 
large time scales, in accord to the two time-scale dynam- 
ics. 

This requires some modification in the simple RIN 
models. The sole use of the delta-correlated Gaussian 
white additive noise, with fluctuating intensity depend- 
ing on u, and a linear force, F(a) = —a, with fluctuating 
7 = 7(1*), is not capable to explain the stretching in the 




-60 -40 -20 20 40 60 



FIG. 7: Conditional probability density function P(a\u) 
given by Eq. for a = e u , k = 1, D = 1.130, B = 0.163, 
Uq — 2.631. The inner curve: u — 0. The outer curve: u = 3. 
x = a/(a 2 ) 1 ^ 2 denotes normalized acceleration. 

observed P(a\u), as it implies only Gaussian conditional 
probability density function P(a\u), for any fixed u. 

However, it is known that accounting for the mul- 
tiplicative delta-correlated Gaussian white noise in the 
drift term of Langevin equation implies stretched expo- 
nential tails. 

Hence we can simply follow the dLDN anzatz as a 
constitutive model (see Sec. IllCf) using the assumption 
that the additive noise intensity a appearing in the sta- 
tionary probability distribution P(a\D, a, B, v$) given by 
Eq. I|65|) depends on u. 

Here, we put the parameter A measuring coupling of 
the noises to each other to zero ignoring thus the skew- 
ness effect which is very small for both the experimental 
P(a\u) and P(a) (l7],|2(j]. This effect is nevertheless of 
much interest since it is associated to the relationship 
between stretching and vorticity in a three-dimensional 
flow. More important here is that it may imply additional 
stretching of the tails as well (see A dependent terms in 
Eq. H A3 ji b This possible way to explain stretched ex- 
ponential tails of the observed P(a\u) can be considered 
elsewhere. 

Following the arguments and techniques presented in 
Sec. Ill A 31 the marginal probability distribution P(a) = 
P(a\D, B, i/q) is obtained by integrating out u in 

P(a\u)=P(a\D,a(u),B,v ), (67) 

with an appropriate choice of the function a(u), for ex- 
ample, a(u) = e u , and some probability distribution of 
u, for example, a Gaussian one with zero mean. Note 
that once P(a\D, cv(it), B, uq) is fitted to the experimental 
curves of the conditional P(a\u) there formally remains 
only one parameter to be fitted in the marginal P(a), the 
variance of the Gaussian distribution of u, i.e., the rms 
velocity u. 

The normalized conditional distribution (|67|l given by 
Eq. f)5J| with a = e u is shown in Fig.0 for four values, 
u = 0,1,2,3, and the other parameters fixed. One can 
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FIG. 8: Conditional probability density function P(a\u) 
given by Eq. JBSJ for a = 1, k = 1, D = 1.13e u , B = 0.163, 
Uo = 2.631. The outer curve: u — 0. The inner curve: u = 3. 
x = a/{a 2 ) 1 ' 2 . 



observe an increase of the variance with the increase of u, 
and a good qualitative agreement with the experimental 
curves P(a\u) 20]. We note that the velocity fluctua- 
tions are present only in the definition of the LDN addi- 
tive noise Q45H which can be viewed as a hint that only 
additive noise intensity essentially depends on u. These 
results give an independent support to the model repre- 
sented in Sec. Ill A 31 

We have checked a different reasonable assumption 
that the multiplicative noise intensity parameter D de- 
pends on the velocity fluctuations, D ~ e", with the 
other parameters fixed. The result is shown in Fig. [S] 
One can observe that the change of the shape of P(a\u) 
defined as P(a\D(u), a, B, vq) with the increase of u 
does not qualitatively meet that observed in the experi- 
ments [20J. We note that an increase of D, i.e., stronger 
multiplicative noise, is generally understood as the pro- 
nounced increase of the relative chance for a fluid particle 
to have higher accelerations as compared to low acceler- 
ations, due to the multiplicative random process. This 
point of view is confirmed by Fig. [S] Also, for a com- 
pleteness in Fig. we represent sample dependencies of 
P(a\D, a, B, Pq) on the parameters B and uq. 

Using the dLDN probability function 
P(a\D, ct(u), B, uq) one can compute {a 2 \u) and 
compare the result with the predicted scaling relation 
(|66ll and the experimental data j2jj. A sample plot 
of (a 2 \u) is shown in Fig. where we have used the 
exponential in the form a = e U//,l ° to provide possibility 
for a fitting to the experimental data (boxes). Although 
with the fitted exponent, uq = 3, the result (triangles) 
exhibits a departure to the experiment, qualitatively 
the model implies a correct behavior of the conditional 
acceleration variance. 

To summarize, the observed stretched exponential 
form of the conditional acceleration probability density 
function P(a\u) can be understood within the framework 
of the dLDN model l|65|l due to the effect of the multi- 
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FIG. 9: Conditional probability density function P(a\u) 
given by Eq. (JHSJ ■ Top panel: a = 1, k = 1, D = 1.130, 
B — 0.163e", vo — 2.631 (the outer curve: u = 0, the inner 
curve: it = 2). Bottom panel: a = 1, k = 1, D = 1.130, 
B — 0.163, vq = 2.631e" (the outer curve: u — 0, the inner 
curve: u = 3). x = a/(a 2 ) 1//2 . 
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FIG. 10: The normalized conditional acceleration vari- 
ance, (a 2 \u) as a function of normalized velocity fluctuations 
u/(u) 1//2 . Boxes: the experimental data on {a 2 \u} / (a 2 ) |2(il |. 
triangles: (a 2 |it)/(a 2 |0) with P(a\u) given by Eq. 165H for 
a = e u/3 , k = 1, D = 1.130, B = 0.163, u = 2.631. 



plicative noise under the assumption that the additive 
noise intensity a depends on velocity fluctuations u. The 
alternative assumption that the multiplicative noise in- 
tensity depends on u seems not to be in a qualitative 
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agreement with the shapes of experimental curves at the 
different values of u except for u = 0. The predicted 
conditional acceleration variance {a 2 \u) with a = e™/ uo , 
uq = 3, have been found in a good qualitative agree- 
ment with the experimental curve. However, we observe 
a departure to the experimental data. 



DISCUSSION 



integral entering the probability density function l|48l) are 
given below. 
At v t = v 0l 



v k 2 a-Da + X D + is Q k 2 2 

da——: = — m\Da — 2Xa + a\ 

Da 2 -2\a + a 2D 1 J 

Xv$k 2 Da — A . . 

arctan (Al) 



ZV-Da - A 2 V-Da - A 2 



We conclude with a few remarks. 

(i) It is interesting to note that a "universal" proba- 
bility density function f(u) could be identified with the 
help of the relation, 



du u e f{u) ~ u 9/2 , 



(68) 



stemming from a comparison of the Heisenberg-Yaglom 
scaling relation l|29|) and the scaling l|t)6|l . both recently 
confirmed by the experiments, where we assume that 
J^^du (a 2 \u)f(u) ~ (a 2 )- Obviously, the choice of a 
normalized Gaussian probability density function with 
zero mean for f(u) does not lead to the above relation 
since it implies ~ u 6 . This can be viewed as a signal for 
selecting a different probability distribution of u, such as 
that of a stretched exponential form. 

(ii) In the present paper, we put the parameter A, 
which measures the coupling between the multiplicative 
and additive noises in the dLDN model, to zero discard- 
ing thus skewness effects in the predicted P(a\u) and 
P{a). This effect is of much interest to study in order 
to estimate A using the fitting to Lagrangian experimen- 
tal data on the longitudinal component of acceleration 
with respect to the trajectory (A = by definition for 
the transverse component of acceleration). Presumably 
it is small due to small skewness of both the observed 
P(a\u) and P(a) for the measured x component of accel- 
eration [l8j . 

(iii) Following the RIN approach presented in this pa- 
per it is of interest to evaluate the proposed averaging 
of the dLDN probability distribution P{a\D, a(u), B, Vq) 
given by Eq. I|65|) over normally distributed u with a 
taken to be a = e u . This is equivalent to the averag- 
ing over log-normally distributed a. A comparison of the 
resulting distribution with the Lagrangian experimental 
data can be made elsewhere. 

(iv) In the presentpaper we have not reviewed a recent 
work by Reynolds A comparison of the results of 
the Reynolds model with that of the stochastic models 
proposed in 0, 01 can be found in the recent paper by 
Mordant et al. 201. 



At vt = B\a\/k, for positive and negative a, respec- 
tively, 



, ^Bka 2 - Da + A Bka 
da — „ — = =p- 



Da 2 -2Xa + a 



D 



D 2 ± 2B\k 
2D 2 



ln[Da 2 - 2Aa + a] 



B(Da-2X 2 )k Da-X 
± , arctan- 



D 2 y/Da — X 2 



y/Da - X 2 



(A2) 



In the general case, we have obtained a cumbersome 
expression, 



da 



—v t k 2 a — Da + X v t k 2 1 , ^ , 



Da 2 -2Xa + a D 2 



-ln[£>a -2Aa + a] 



2BXk 
~D 2 ~ 



\n[2Bka + v t k 2 ] + F{c) + F(-c), (A3) 



where we have denoted 



2c 2 D 2 c cic 2 (c-Da + X) 
x (B 2 (A 2 + cA - Da)a + c(Dv 2 k 2 + c 2 v t ))], (A4) 



c = —iy/ Da — A 2 , v t = Jvl+ B 2 a 2 /k 2 , (A5) 



Cl = S 2 (4A 3 + 4cA 2 - 3-DoA - cDa) + D 2 (c + X)is$k 2 , 

(A6) 



c 2 = \ B 2 (2X 2 + 2cX - Da)k 2 + D 2 v 2 k A . (A7) 



Some useful formulas used in verifying the limits B — > 
and D — > are: 



arctana; = ^(^(-^ ~ — m (l + * a; ))' (^) 



APPENDIX A: EXACT INTEGRALS 

Exact indefinite integrals, up to a constant term which 
does not depend on a, used in calculating the definite 



lim — lnfl + Da 2 } = a 2 . 
d^o D 



(A9) 
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